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[ — . , Abstract 

We consider the fluid interface problem given by two incompressible fluids with dif- 
ferent densities evolving by Darcy's law. This scenario is known as the Muskat problem 

■ for fluids with the same viscosities, being in two dimensions mathematically analogous to 
the two-phase Hele-Shaw cell. We prove in the stable case (the denser fluid is below) a 

■ maximum principle for the L°° norm of the free boundary. 

1 Introduction 



The Muskat problem models the fluid interface problem given by two fluids in a porous 
medium with different characteristics. The problem was proposed by Muskat (see [13J) in 
a study about the encroachment of water into oil in a porous medium. In this phenomena, 
Darcy's law is used to govern the dynamics of the different fluids [2j. This law is given by 
£NJ ■ the following formula: 



o 



-v = -Vp-(0,0,gp), 

K 



where v is the velocity of the fluid, p is the pressure, p is the dynamic viscosity, k is the 
permeability of the isotropic medium, p is the liquid density and g is the acceleration due to 
^ ■ gravity. The fluid interface is given between incompressible flows with different viscosities, 

pi, fj,%, and densities, pi, p%. 

Saffman and Taylor [14J considered this problem in a study of the dynamics of the interface 
between two fluids with different viscosities and densities in a Hele-Shaw cell. In this physical 
scenario (see [11]) the fluid is trapped between two fixed parallel plates, that are close enough 
together, so that the fluid essentially only moves in two directions. The mean velocity is 
described by 

_^ = _Vp-(0,gp), 

where b is the distance between the plates. Darcy's law, in two dimensions, and the above 
formula become analogous if we consider the permeability of the medium k equal to the 
constant b 2 /12. 

The Muskat problem and the two-phase Hele-Shaw flow have been extensively considered 
(see [1] and [12] and the references therein). These free boundary problems can be modeled 
with surface tension [9] using the Laplace- Young condition. In this case there is a jump 
of discontinuity in the pressure of the fluids across the interface proportional to the local 
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curvature of the free boundary. In this work we will consider the case without surface tension, 
so that the pressures are equal on the interface. 

This scenario is considered by Siegel, Cafhsch and Howison in [15] . where they show ill— 
posedness in an unstable case and global-in-time existence of small initial data in a stable 
case. They describe the two-dimensional dynamics of the incompressible flow as follows: 

v = -aVp-(0,V), 

where a takes two positive constant values 

for i = l,2, 



12 W ' 

on each fluid, and V is a constant. With our notation, this case is equivalent to consider in 
the two-dimensional problem 

mi, v 

— = — , and g Pi = —, 

Ki (Xi (X% 

for % = 1,2. The results rely on the assumption that the Atwood number 

Ml - M2 



Ml + M2 



is nonzero, and therefore there is a jump of viscosities. 

In the same year, Ambrose [T] studies the 2-D problem with initial data fulfilling 

(pa - Pi)g cos(6(a, 0)) + 2A lx U(a, 0) > 0, 

where 9 is the angle that the tangent to the curve forms with the horizontal and U is the 
normal velocity (given by the Birkhoff-Rott integral). In this work he considers the arclength 
and the tangent angle formulation used by Hou, Lowengrub and Shelley [12J to get energy 
estimates for the free boundary assuming the geometric condition 

(x(a,t) -x(a',t)) 2 + (y{a,t) -y(a',t)) 2 
[a — a'Y 

locally in time, where the curve (x(a,t),y(a,t)) is the interface. We show in |10j that this is 
not enough to get local-existence in this kind of contour dynamics equations, since a regular 
interface could touch itself with order infinity and without satisfying the inequality ([1]) . 

We will study the fluid interface given by a jump of densities, so that = 0. Therefore 
we can take k = b 2 /12 = 1, and fx\ = ^ = g = 1 without loss of generality. This case 
can consider, among others, the dynamics of moist and dry regions in porous media. This 
scenario was considered by Dombre, Pumir and Siggia [8j, but in a different context. They 
treated the interface dynamics for convection in porous media where, using our notation, the 
density plays the roll of the temperature. They studied the unstable case, namely when the 
denser fluid (or the fluid with larger temperature) is above. They present meromorphic initial 
conditions with complex poles, and study the dynamics of these critical points. 
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It is well-known [12J that for these contour dynamics systems, the velocity in the tangential 
direction does not alter the shape of the interface. If we change the tangential component of 
the velocity, we only change the parametrization. In [7J we used this property to parameterize 
the interface as a function (x, f(x,t)), getting the following equations: 

f t (x t) = P^PV [ ^f(^)-Vf(x- y) t)).y 

J{ ' ' 4tt i K2 [\y\* + (/(*, t) - f(x - y, t)Yf/z y ' (2) 

f(x,0) =/ (x), xeK 2 , 

for a two-dimensional interface, and 

f( +\- P2 ~ pl pv f (j^/fa l l ~ dx f( a ~ /3, * ^ ^ 

M a,tj- ^ ^y Rj g a + ( / ( 0|t) _ / ( _ ) g >t )) 2 ^ (3) 

/(a,0) = /o(a), ael, 

for a one-dimensional interface. We point out that with these formulations the condition (pTJ) 
is satisfied locally in time if local-existence for the systems is reached. This avoids a kind 
of singularity in the fluid when the interface collapses (see [6] for example). Also we prove 
that when the denser fluid is below the other fluid, p 2 > Pi, the problem is well-posed given 
local-existence and uniqueness for the systems ([2]) and ([3]). When the less dense fluid is below, 
pi < pi, we prove ill-posedness showing that the equations ([2]) and © are ill-posed. We get 
this result using global solutions of ((3j) in the stable case, p% > pi, for small initial data in a 
similar way as in [15] . 

If we neglect the terms of order two in ([2]), the linearized equation is obtained. It reads 

ft = ^^(Rid Xl f + R 2 d x J) = ^-^A/, 
f(x,0) = f (x), 

where R\ and R 2 are the Riesz transforms (see [16] ) and the operator A/ is given by the 
Fourier transform A/(£) = In the stable case pi < p 2 (the greater density is below), 

the linear equation is dissipative and is clear that the following maximum principle is reached 

[|/[|r»(t) < [|/o[|r». 

In this work we get, in section 3 and 4, this estimate for the nonlinear systems (j2J and ([3]). 
To this end, we follow the evolution of the maximum of the absolute value of f(x,t). This 
technique was used for one of the authors in [5] in a family of dissipative transport equations 
for incompressible fluids. Also we would like to cite the work of A. Constantin and J. Echer 
where they study the shallow water equation in the same way. By a similar approach, in 
section 5 we obtain a global bound on the derivative for small initial data. 



2 Parameterizing the interface in terms of a function 

In this section we briefly explain how to parameterize the free boundary in terms of a function 
(see [7] for more details). This way of writing this nonlocal equation is crucial to check the 
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evolution of the maximum of the absolute value of the function to obtain the maximum 
principle. 

In our case, Darcy's law can be written as follow: 

v(xi,x 2 ,x 3 ,t) = -Vp(xi,x 2 ,x 3 ,t) - (0,0,/j(xi,x 2 ,x 3 ,t)), (5) 

where {x%, x 2 , x 3 ) £ 1R 3 are the spatial variables and t > denotes the time. Here p is defined 
by 

pi in Qi(t) 



p(x U X2,X 3 ,t)- l ^ . n ^ 

with pi,p2 > constants and pi ^ pi. The set is given by 

n x (t) = {x 3 > f(xi,x 2 ,t)}, 

and 

n 2 (t) = {x 3 < f(xi,x 2 ,t)}, 

being f(x\,x 2 ,t) the fluid interface. If we apply the curl operator to Darcy's law twice the 
pressure disappears. Considering the incompressibility of the fluid, we have curl curl v = — Av, 
and we can give the velocity in terms of the density as follows: 

v = {d Xx /ST 1 d X3 p,d X2 A~ 1 d X3 p : -d xl A~ l d xl p - d X2 A~ l d X2 p). (6) 

The density p has a jump of discontinuity on the free boundary, therefore the gradient of it 
function is given by a Dirac distribution 5 as follows: 

Vp= (p 2 - p 1 )(d Xl f(xi,x 2 ,t),d X2 f(xi,x 2 ,t),-l)5(x 3 - f(xi,x 2 ,t)), (7) 

Using the kernels for d Xl A _1 and d X2 A~ l we obtain 

/ ,x p 2 -px [ (2/1,2/2, V/(x - y,t) ■ y) 
V (xi,X 2 ,X3,t) = T— PV n | 2 77 ^213/2 ^' 8 

where we note x = (xx,x 2 ), y = (yi,y 2 ). In ([8]) x% ^ f(x,t) and the principal value is taken 
at infinity (see [E]). The vorticity is at the same level than the gradient of the density, so 
it is given by a delta function. This forces the velocity to have a discontinuity on the free 
boundary. Just checking the incompressibility of the fluid in the sense of the distributions, 
we obtain that this discontinuity is in the tangential direction, so that it does not affect the 
shape of the interface (see [7])- Ignoring the tangential terms we obtain that the velocity on 
the free boundary is given by 

v(x fix t) t) = - p2 - pl py f iyuy2,Vfjx-y,t).y) 

( ' j 4vr FV L[\y\i + (f(x,t)-f(x-y,tm^ dy - (9) 

If we want to parameterice the free boundary in terms of a function, we need to get the 
velocity v = (i>i,t>2> ^3) satisfying v\ = v 2 = 0, since otherwise the points on the plane are 
not fixed and they depend on time. If we add the following tangential terms to ([9]): 

02 ~ Pl - PV I Fi ,9 , T, V \, —^dyil,0,d Xl f(x,t)), 



4vr 7 R2 [\y\ 2 + {f{x,t)-f{x-y,t)Yf/ 2 
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El T^ PV n i2T7T( \ -^dy(0,l,d x J(x,t)), 

4vr 7 R 2 + (f(x, t) - f(x - y, t)) 2 }^ 2 

we do not alter the interface and this follows: 

„(„, /<,, *),*) = «_« to, o, pv f gf^^pML^ ,,). (10 , 

Finally the contour equation given by 

f( t )- P2 - Pl py [ (Vf(x,t)-Vf(x-y,t))-y 

Jt{ ' ' " 4vr i R2 [|y|2 + (/( x , t ) - f(x - y, t)) 2 f/ 2 V ' (11) 

/(x,0)=/ (x). 

This formula works for periodic interface and for a free boundary near planar at infinity. In 
both cases it presents an principal value only at infinity. If we suppose that the function 
f(x,t) only depends on x\, integrating in a;2, the contour equation in the 2-D case is reached 

f(«.*\- p2 ~ pl PV f (d x f{x,t)-d x f(x-a,t))a 

M x,t)- ^ rv J na 2 + {f{X:t) _ f{x _ ajt)) 2 aa > (12) 

f(x,0) =/ (x); xeR. 

This equation can be obtained in a similar way that (lllh using the stream function. We check 
in [7J that as long as this equation is satisfied we obtain weak solutions of the system 

p t + v ■ Vp = 0, 

(13) 

v = -Vp- (0,0, p), divw = 0. 



3 Two dimensional case (1-D interface) 

Here we show that the L°° norm of the system (|12p decreases in time in the stable case 
(p2 > pi). We will consider the set O equal to R or T. The following theorem is the main 
result of the section. 

Theorem 3.1 Let /o G H k (Q) with k > 3 and p2 > p\. Then the unique solution to the 
system l\12\\ satisfies the following inequality: 

\\fh~(t) < ||/olU«. 

Proof: For /q £ H k with k > 3, we prove in [7J that there exists a time T > so that 
the unique solution to (fT2j) f(x,t) belongs to C 1 ([0, T]; H k ). In particular we have f(x,t) S 
C 1 ([0,T] x 0), so that the Rademacher theorem gives the functions 

M(t) = max/(x,t), 

X 

and 

m(t) = mm f(x,t), 
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differentiable almost every t. In the non periodic case, we notice that also there always exists 
a point ij G K where 

\f(x t ,t)\ = max\ f(x,t)\, 

X 

due to the fact that f(-,t) 6 H s with s > 1/2, and using the Riemann-Lebesge lemma, f(x, t) 
tends to when \x\ — > oo. We suppose that this point x% satisfies that < f(xt,t) = M(t). 
A similar argument can be used for m(t) = f(xt,t) < 0. If we consider a point in which M(t) 
is differentiable, we have 

= Km M(t + h)-M(t) 
h-^o+ h 

f{x t +h,t + h) - f{x t ,t) 
h->0+ h 

,. f(x t+h ,t + h)-f(x t ,t + h) f(x t ,t + h) - f(x t ,t) 

= hm 1 . 

h->o+ h h 

Since f(x, t + h) takes on its maximum value at x = x t+ h, it follows: 

wz/.w ,. f{x t ,t + h) - f(x t ,t) 

M (t)> hm = f t (x t ,t). 

h-*o+ h 

Computing for /i > 

M(t)-M(t-h) 

M'(t) = lim — — i - 

h^o+ h 

_ y f{x t ,t) - f{x t -h,t- h) 

h->o+ h 

y m f(xt,t- h) - f(x t -h,t- h) f(x t ,t) - f(x t ,t- h) 
7i->0+ h h 

< lim f{xt,t) - f{x t ,t- h) 
~ h-*0+ h 

< ft{x t ,t), 

and we obtain finally 

M'(t) = f t (x t ,t). (14) 
If we take the value x = xt in the equation (|12h . the above identity gives 

W 2vr J RO ? + {{f{x t ,t)- f{x t -a,t))f 

using the fact that d x f(xt,t) = 0. Integrating by parts 

M'(t) = P2 - Pl PV [ ^(/(st.*) ~ f^ Xt ~ I da 

2vr J R a 1 | ^ f(xt,t) - f(x t - a,t) y a 

= h + h, 
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where 



j _ P2 ~ Pi p y f f(x t ,t) - f(x t - a,t) . ^ 

2vr J R a 2 1 ff(xt,t)-f(xt-a,t)\2 



a 
and 

( f(x t ,t) - f(x t - a,t) y 

h = -iLi£L [ 2 - \ fl ° i— sj l^ zfia -"■') ) 



a 



Using the function 

we can write I2 as follows: 



Gfx) = K + arctanx, 

1 + x z 



I 2 = -^^PV [ d a G( f{Xut) ~ nXt ~ a ^)da. 
2vr J R \ a ) 

Integrating we obtain 
The /1 term is equal to 

2ir y R a 2 + (M(t)-/(x ( -a,t)) 2 - 
so that M'(t) < for almost every i. In a similar way we obtain for m(t) the following: 

m'{t) = -E^PV [ m(t)- f( Xt -a,t) 

W 2vr 7 M a 2 + (m(i) - /(^ -a,i)) 2 " 

for almost every t. Integrating in time we conclude the argument obtaining the maximum 
principle. 

Let = T, with this maximum principle we can conclude the following decay of the L°° 
norm. 

Proposition 3.2 Let /o G H k (T) with k > 3 and pi> P\- If 

fo{x)dx = 0, 



T 

then the unique solution to the system (|12p satisfies the following inequality: 

||/IU-(*)<||/o|U-e- ( ' ,2 -' ,l)C ' (ll/o|l£0c) *, 

with C(||/ ||ioc) > 0. 



Proof: Suppose that 

( f (x)dx = 0. 

We can write (|12|) as follows: 



2tt I k V a 



and therefore we have 



f t {x,t)dx = f PV [ ^.arctanf ^-^-^W x 

_P2~ Pi pv f f d , an ( f(x,t)-f(x-g,ty 



' 2 ~ JrJi v - « 
= 0. 

Integrating in time we obtain 

/ f{x,t)dx = 0, Vt > 0. (15) 

As we show in the proof of the above theorem, we have 

d imi n\ ft-Pi mrf \\fh~>(t) - f(x t - a,t) 

dt mL °°v = -— pv J Ra l +mLOO{t) - fixt - a , tw da > 

for almost every t. Using the maximum principle, for \a\ < r we get 

« 2 + (\\fh~(t) - f(x t - a,t)) 2 <r 2 + 411/ollicc, 

and it follows: 

d imi u\*> P2 ~ Pi mr f \\f\\L°°(t)- f(x t -a,t) 

at 2vr J\a\<r a +{\\f\\L°°{t)- f{x t - a,t)Y 

p 2 -pi 2r 92- P\ 1 (tt \j 

o 2 i /in f 112 + — 5 2 , ,11, i|2 / f(x t -a)da. 

2vr H + 4 /0 2,00 2?r r + 4 /o !/*> /|a|<r 



If we take r = n-zr for n G N, using (|15p we obtain 

£n/iu-(t) < - V 1 2 2 +"I,,, 2 ii/iu-c). 

and integrating in time we conclude the proof. 
For $7 = K we obtain the following result. 

Proposition 3.3 Let /o G mf/i > 3 and P2> p\- If fo(%) < or /o(x) > ; then 

the unique solution to the system (|12p satisfies the following inequality: 



l + (p2-pi)C(||/ ||^,||/o|| L i)r 

with C(||/o||ioo,||/o|| z i) > 0. 



8 



Proof: Let consider fo(x) > 0. The argument is similar in the other case. Our maximum 
principle shows that 

,,,, P2-P1 /" m{t) - f(x t - a,t) 

m (t) = PV / — = ; — ^ —. TT^da > 0, 

W 2vr 7m a 2 + (m(t) - /(z t - a, t)f ~ 

for almost every t, so that if fo(x) > 0, then f(x,t) > 0. In a similar way as in the previous 
result, we can conclude that 

/ f t (x,t)dx = 0, 

and therefore 

/ f(x,t)dx = / f (x)dx. 

Due to / is nonnegative, we control the L 1 norm of the solution, so that ||/||x,i(t) = H/oIIl 1 - 
We have = f(xt,t), and 

^II/I|l-(*) = -/, 

with 

/ = Pl^flpv I f(xt ; t)- f(x t -a t t) 



2vr J R a2 + (f(x t ,t)- f(x t -a,t)) 2 

for almost every t. If we consider the interval [— r, r] for r > 0, 

tfi = {« G [-r, r] : /(x t , t) - f(x t - a,t) > f(x t , t)/2}, 

and 

U 2 = {a € [-r, r] : /(x t , i) - /(x t - a, i) < /(x t , t)/2}, 

we get 

T> P2-P1 pv f f(x t ,t) - f(x t - a,t) tJn> P*Zl£l f^t)/ 2 \tj I 

" 27T 7 [/ia 2 + (/(xt)t) _ /(xt _ a)t)) 2 - 2?r r 2 +4 ||/ ||2 oo l X| - 

In order to estimate \Ui\, we use that \U\\ = 2r — jl^l, and 

/• /• f(xt,t) 

II/oIIli = / f{x t - a,t)da> / f(x t - a,t)da> — |J7 2 |, 

and therefore > 2(r — ||/o||li *))• This estimate gives 

T P2-P\ f(x t ,t)/2 92 - P\ rf(x t ,t) - \\f \\ L i 

~ 2tt r 2 + 4 || /0 ||| oo |tyi1 - 2tt r2 + 4||/ ||i« ' 

and this function reaches its maximum at 



r = (II/oIIli + J\\M\ 2 L1 +4||/ ||| oo / 2 (x t ,t))//(^,t). 
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Using the maximum principle, easily we get 



P i ~ Pi \\fo\\L^f 2 (xt,t) 

8tt ||/o||! 1 +2||/o|Ui||/o|l!oc+2||/ ||| a 



1 > "o-'h* ,12 i 1 J <TL K 7^\ 7^1 „4 > (P2 " Pl)C(||/o|kx, \\h\\L~)f 2 {xut). 



Finally, we obtain 

d 



(// „., ^At) < -(P2 - pi)c(||/o||xi, [|/olk«)[|/[|i»(t), 



and integrating we end the proof. 



4 Three dimensional case (2-D interface) 

In this section, by using the same technique, we extend the maximum principle for the three 
dimensional stable case. We consider the set f2 the plane or the periodic setting. 

Theorem 4.1 Let /q G H k (Q) for k > 4, and p2 > pi- Then the unique solution to (fTTj) 
satisfies that 

< ||/o||l°°- 



Proof: As we prove in [7J, there exists a time T > and a unique solution f(x,t) G 
C 1 {[0,T};H k (n)) solution of PJ. In particular f(x,t) G C^QO.T] x n) using Sobolev in- 
equalities. In the case Q = M 2 , there always exists a point xt G M 2 where \ f(x,t)\ reaches 
its maximum due to the fact that f(-,t) G H s with s > 1, and using the Riemann-Lebesgue 
lemma f(x, t) tends to when |x| —> oo. Suppose that this point is for M(t) = f(xt, t) > 0. A 
similar argument can be used for m{t) = f(xt,t) < 0. By using the H. Rademacher theorem, 
the function M(t) is differentiable almost everywhere, and computing as before we obtain 

M'(t) = f t (x t ,t), (16) 

for almost every t. Using equation (fTTj) . the fact that Vf(x t ,t) = 0, and the last identity, we 
have 

M '( f \ = P2-Pi pv f -Vf(y,t) ■ (x t -y) 

{ ' 4tt J R2 [\ Xt - y| 2 + (f(x t , t) - f(y, t))^ V ' 

Integrating by parts 

M\t) = El^lpvf VyU{xut) _ /(y>t)) . JiLiyfr + f f(*ut)-f(v,t) yyw 

4vr J R 2 \x t -y\ 6 \ V \x t -y\ J J 

= -HLZHLrvf {fM - fM)(iiVy J!LZj' ) l 1 + ( /(»■■«) -/(».0 )V'^ 
4vr 7 R2 kt-y^'V V |a; t - y\ J J 

P2-Pi p yf f(x t ,t) - f(y,t) x t -y v / | / /(g f ,t) -/(j/, t) \ 2 \- 3 / 2 d 

4vr 7 R 2 kt-yl kt-y| 2 ^ \ \x t -y\ > ' 
= Ji + J 2 . 
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We have 

t P2-Pl DT/ /"«/,, in rrf f(Xi,t) ~ f(y,t) 

J 2 = / V^ln x t - y ) • V y H[ , , 



where 



;r 3 



1+x 2)3/2- 

The identity A 3/ (ln|xt — y\)/4m = 6(xf), and the following limit: 

Hm f(x t ,t)-f(y,t) = Hm nxt,t)-f(y,t)-Vf(x t ,t)-(x t -y) = Q 
y^x t \xt — y\ y-^x t \xt — y\ 

show that integrating by parts in J 2 , we obtain 

J 2 = P -^PV [ A y (ln\ Xt - y\)H( ^f ~ f ^ )dy = ( P2 - Pl )H(0), 

and therefore J 2 = 0. The Ji term is equal to 

P2-P1 f M(t)-f(y,t) 
1_ 4tt y Ra [|x t -i/| 2 + (M(t)-/(y,OT/! 1 y ~ ' 

so that M'(t) < for almost every t. For m(i) we have m'(t) > 0. Integrating in time we 
conclude the proof. 

As in the previous section, using this maximum principle we get the following decay of 
the L°° norm. 

Proposition 4.2 Let f G H k (T 2 ) with k > 4 and p 2 > pi- If 

fo(x)dx = 0, 



It 2 

then the unique solution to the system (|12p satisfies the following inequality: 

||/||^-(t) < H/olU-e-C^-pOCCII/QlU-)*, 

with C(||/o|Uoc) > 

Proof: We can write (jlip as follows: 



47T y K 2 \y\^ \ \y> 



f(x,0) = f (x), 

with 



P(x) 



y/l + X 2 
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As in the previous section, checking the evolution of the integral of / on T 2 , we obtain 

/ f(x,t)dx = 0. (17) 

The proof in the above theorem shows that 

j*ll f H m _ P2 ~ Pi f f ll/IU°°(*) - f(y,t) 

dt MU ™ [t) ~ 4vr 7 R2 [\x t - y\ 2 + (\\f\\L~(t) - /M)) 2 F 2 ^ 

for almost every t. If we consider xt — y G [— rwr, mr] x [— n-/r, n7r] = j4„, with n G N, we have 

\x t - y\ 2 + (\\fh<~(t) - f(x t - a,t)) 2 < 2(nvr) 2 + 4[|/ ||ioo. 
Using (fT7]) . the above estimate gives 

d ii#n m^- P2 - pi mr f ll/IU°°(*) - f(y,t) , 

dt mLA t] ~ ~^r PV J {xt - y) eA n [\ Xt -y\* + mLao{ t)-ny,t))*]W dy 
P2-P1 (2mr) 2 
* ~ 4. [2^^1«^F II/Ik00( ) - 
Integrating in time we finish the proof. 

Proposition 4.3 Let f G H k (R 2 ) with k > 4 and p 2 > p\. If f (x) < or f (x) > 0, then 
the unique solution to the system (|12p satisfies the following inequality: 



^ \\Jo\\L= 



(l + (p2-pi)C(||/ ||L-,||/o|ki)t) 2 ' 

with C(||/o|U-, H/olUO > o. 

Proof: Let us consider fo(x) > 0. The same estimate is obtained for fo(x) < 0. We obtain as 
before f(x,t) > 0, and ||/|| x i(t) = ||/ || L i. We have \\f\\ L ™(t) = f(x t ,t), and 

J t \\fh°°(t) = -J, 

for almost every t, with 

T _ Pi- Pi p y f fjxut) - f(y,t) 

^ Jwz[\xt-y\ 2 + (f(x t ,t)-f(y,t))*}yz y - 

If we define the set B r (xt) = {y ■ \xt — y\ < r} for r > 0, 

Vi = {y G B r (s t ) : /(a; t ,t) - /(y,i) > /(z t ,t)/2}, 

and 

^2 = {y G fl r (a*) : - f(y,t) < f(x t ,t)/2}, 
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we get 

j > P2-P1 f(xt,t)/2 
- 4vr [ r 2 + 4 || /o ||2 oo] 3/2l H 



Using that | V\ | = nr 2 — \ V2 | , and 



||/o||xi> / f(y,t)dy>^^\v 2 \, 
Jv 2 1 



this estimate follows: |Vi| > ixr — 2\\fo\\ L i / f(xt, t). We have 

T > 92- P\ nr 2 f(xt,t) -2||/ || L i 

" 8tt [ r 2 +4 || /o ||2 oc] 3/ 2 

and taking 



2||/oIIli/vt + 1\1/2 



f(xt,t) 
we find 

j > P2ZPI <f{Xi,t)f/ 2 > P2 - Pl (f(x t ,t)f/ 2 



8^ [l + 2||/o[UiA + 4||/ ||i<»/(ar t ,*)]3/3 " 8 [1 + 2||/ |Ui/7T + 4||/ ||3 oo ]3/2 • 
Finally, the following estimate is obtained: 

J t \\fh-(t)<-(p2-pi)C(\\f \\ L i, \\fo\\L~)\\ff£(t), 

and integrating we end the proof. 

5 Small initial data 

In the two-dimensional case, we prove in [7] that if the following quantity of the initial data 
is small 



then there is global- in-time solution of the system (|12p . Here we show that if initially the 
L°° norm of the first derivative is less than one, it continues less than one for all time. 

Lemma 5.1 Let fo £ H s with s > 3, and \\d x fo\\L°° < 1- Then the unique solution of the 
system (|12|) satisfies 

\\d x fh°°(t) < 1. 



Proof: If we consider the following term in (| 1 2 j) : 



2ir J R a 2 + (f(x,t)-f(x-a,t)) 2 
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we can integrate by parts and get 

R _ P2 ~ pi p y r d a (f(x,t) - f(x- a,t)) 1 ^ 

2vr 7 R a 1 / f(x,t) - f( x - a,t)\ 2 

V a 

= -p^fipv f /(*.*)-/(*-<».*) — 1 da 

2vr J R a 2 1+ (j{x,t)-}{x-a,t)\ 2 



a 

' f(x,t)-f(x-a,t) y 
92- Pi f n [ a / q ( f(x,t)- f{x-a,t) ,^ 



a 

Li+L 2 . 



As we showed before 



- f(x-a,t) y 



L 2 = - P2 Pl PV / d n G( Jy ~'"' JV " ~'" / )da = 0. 



2vr i R V a 

so K = L 2 . Making a change of variables we find the following equivalent system: 

,/ ,s P2 ~ Pi m/ f d x f(x, t)(x -a)- (/(a, t) - f(a, t)) 

ft{x,t) = — PV / da. 

2vr 7 R (x - a) 2 + (f{x, t) - /(a, t)) 1 

Taking one derivative in this formula, we have 

d x f t (x)=N 1 (x)+N 2 (x), (18) 



with 



where 



and 



2vr y R (x-a) 2 + (/(a;)-/(Q)) i 
2vr 7 R (x - a) z 



A a /(x) 



(l + (A a /(x))2)5 

/(x) - /(a) 



x — a 

Now, we define in this section 

M(t) = \\d x f\\ L ~{t), 

then M(i) = max x d x f(x, t) = d x f(xt, t) where xt is the trajectory of the maximum. Similar 
conclusions are obtained for m(t) = mm x d x f(x,t). Using the Rademacher theorem as in 
the previous section, we have that M'(t) = d x ft(xt,t), and d%f(xt,t) = 0. Therefore taking 
x = xt in (JIB]), we get 

M'(t) = N 2 (x t ), 
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due to Ni(xt) = 0. The inequality 



|A a /(x t )| < M(t), 

shows that for M(t) < 1 the integral N 2 (x t ) < 0, and therefore M'(t) < 0. If M(0) < 1, 
using the theorem of local existence, we have that for short time M(t) < 1, and therefore 
M'(t) < for almost every time. This implies that M{t) < 1. For m(t) we find m(i) > 1. 
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